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This paper aims at advancing the field of electro-sensing. It exhibits 
the physical mechanism underlying shape perception for weakly elec- 
tric fish. These fish orient themselves at night in complete darkness 
by employing their active electrolocation system. They generate a 
stable, high-frequency, weak electric field and perceive the transder- 
mal potential modulations caused by a nearby target with different 
admittivity than the surrounding water. In this paper, we explain 
how weakly electric fish might identify and classify a target, know- 
ing by advance that the latter belongs to a certain collection of 
shapes. The fish is able to learn how to identify certain targets and 
discriminate them from all other targets. Our model of the weakly 
electric fish relies on differential imaging, i.e., by forming an image 
from the perturbations of the field due to targets, and physics-based 
classification. The electric fish would first locate the target using 
a specific location search algorithm. Then it could extract, from 
the perturbations of the electric field, generalized (or high-order) 
polarization tensors of the target. Computing, from the extracted 
features, invariants under rigid motions and scaling yields shape de- 
scriptors. The weakly electric fish might classify a target by com- 
paring its invariants with those of a set of learned shapes. On the 
other hand, when measurements are taken at multiple frequencies, 
the fish might exploit the shifts and use the spectral content of the 
generalized polarization tensors to dramatically improve the stability 
with respect to measurement noise of the classification procedure 
in electro-sensing. Surprisingly, it turns out that the first-order po- 
larization tensor at multiple frequencies could be enough for the 
purpose of classification. A procedure to eliminate the background 
field in the case where the permittivity of the surrounding medium 
can be neglected, and hence improve further the stability of the 
classification process, is also discussed. 

weakly electric fish | electrolocation | shape classification | generalized polar- 
ization tensors | shape descriptors | spectral induced polarization 

Introduction 

In the turbid rivers of Africa and South America, some species 
of fish generate a stable, high frequency, weak electric field 
(0.1-10 kHz, < 100 mV/cm) which is not enough for defense 
purpose. In 1958, Lissmann and Machin [22] discovered that 
the emitted electrical signal is in fact used for active electro- 
sensing. The weakly electric fish have thousands of receptor 
organs at the surface of their skin. A nearby target with dif- 
ferent admittivity than the surrounding water is detected by 
measurements of the electric organ discharge modulations at 
the receptor organs [131 125| . Targets with large permittivity 
cause appreciable phase shifts, which can be measured by re- 
ceptors called T-type units [2B]. It is an important input for 
the fish, and thus it will be the central point in this paper for 
shape classification. 

Active electro-sensing has driven an increasing number of 
experimental, behavioral, biological, and computational stud- 
ies since Lissmann and Machin's work [TT1 ri2l [T5l [TTl [23l [24l 
1301 136] . Behavioral experiments have shown that weakly elec- 
tric fish are able to locate a target [36] and discriminate be- 
tween targets with different shapes [35] or/and electric pa- 
rameters (conductivity and permittivity) [33] . The growing 
interest in electro-sensing could be explained not only by the 
curiosity of discovering a "sixth sense", the electric percep- 



tion, that is not accessible by our own senses, but also by po- 
tential bio-inspired applications in underwater robotics. It is 
challenging to equip robots with electric perception and pro- 
vide them, by mimicking weakly electric fish, with imaging 
and classification capabilities in dark or turbid environments 

[M DS ED L2Z1 LH 131 133] • 

Mathematically speaking, the problem is to locate the tar- 
get and identify its shape and material parameters given the 
current distribution over the skin. Due to the fundamental 
ill-posedness character of this imaging problem, it is very in- 
triguing to see how much information weakly electric fish are 
able to recover. The electric field due to the target is a com- 
plicated highly nonlinear function of its shape, admittivity, 
and distance from the fish. Thus, understanding analytically 
this electric sense is likely to give us insights in this regard 
[IT] EH IH EGD HI H3 [33- While locating targets from the 
electric field perturbations induced on the skin of the fish is 
now understood (see [JJI2JJ and references therein), identify- 
ing and classifying their shapes is considered to be one of the 
most challenging problems in electro-sensing. Although the 
neuroethology of these fish has significantly been advanced 
last years (see |16] and references therein), the neural mecha- 
nisms encoding the shape of a target is far beyond the scope of 
our study. Rather, this work focuses on the physical feasibility 
of such a process, which was not explained until now. 

In p], a rigorous model for the electro-location of a target 
around the fish was derived. Using the fact that the elec- 
tric current produced by the electric organ is time-harmonic 
with a known fundamental frequency, a space-frequency lo- 
cation search algorithm was introduced. Its robustness with 
respect to measurement noise and its sensitivity with respect 
to the number of frequencies, the number of sensors, and the 
distance to the target were illustrated. In the case of disk- 
and ellipse-shaped targets, the conductivity, the permittivity, 
and the size of the targets can be reconstructed separately 
from multifrequency measurements. Such measurements have 
been used successfully in trans-admittance scanners of breast 
tumors [T01I201I3T] . 

In this paper, we tackle the challenging problem of shape 
recognition and classification. In order to explain how the 
shape information is encoded in measured data, we first derive 
a multipolar expansion for the perturbations of the electric 
field induced by a nearby target in terms of the characteris- 
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tic size of the target. Our asymptotic expansion generalizes 
Rasnow's formula [25] in two directions: (i) it is a higher- 
order approximation of the effect of a nearby target and it is 
valid for an arbitrary shape and admittivity contrast and (ii) 
it takes also into account the fish's body. As it has been first 
proved in [T] , by postprocessing the measured data using layer 
potentials associated only to the fish's body, one can reduce 
the multipolar formula to the one in free space, i.e., without 
the fish. Then we show how to identify and classify a target, 
knowing by advance that the latter belongs to a dictionary 
of pre-computed shapes. The shapes considered in this paper 
have been experimentally tested and results reported in [34] , 
This idea comes naturally in mind when modeling behavioral 
experiments such as in |33ll35l[5rJ] . The pre-computed shapes 
would then be a model for the fish's memory (trained to rec- 
ognize specific shapes), and the experience of recognition pre- 
sented here would simulate the discrimination exercises that 
are then imposed to them. We develop two algorithms for 
shape classification: one based on shape descriptors while the 
second is based on spectral induced polarizations. We first ex- 
tract, from the data, generalized (or high-order) polarization 
tensors of the target (GPTs) [2]. These tensors, first intro- 
duced in [6], are intrinsic geometric quantities and constitute 
the right class of features to represent the target shapes 19] . 
The shape features are encoded in the high-order polarization 
tensors. The extraction of the GPTs can be achieved by a 
least-squares method. The noise level in the reconstructed 
generalized polarization tensors depends on the angle of view. 
Larger is the angle of view, more stable is the reconstruc- 
tion, /i-regularization techniques could be used. Then we 
compute from the extracted features invariants under rigid 
motions and scaling. Comparing these invariants with those 
in a dictionary of pre-computed shapes, we successfully clas- 
sify the target. Since the measurements are taken at multiple 
frequencies, we make use of the spectral content of the gen- 
eralized polarization tensor in order to dramatically improve 
the stability with respect to measurement noise of the physics- 
based classification procedure. In fact, we show numerically 
that the first-order polarization tensor at multiple frequencies 
is enough for the purpose of classification. 



Feature extraction from induced current measurements 
Electro-sensing model. Let us recall the nondimensionalized 
model of electro-sensing: the body of the fish is f2 (of size of 
order 1), an open bounded set in M. d ,d = 2,3, of class C 1,Q , 
< a < 1, with outward normal unit vector v. The electric 
organ is a dipole / placed at zq G Q or a sum of point sources 
inside Q satisfying the neutrality condition. We refer to [T] 
where the equations are nondimensionalized and the different 
scales are identified. The fish's skin is very thin and highly 
resistive. Its effective thickness, that is, the skin thickness 
times the contrast between the water and the skin conductiv- 
ities, is denoted by £ and is of order 10 -1 [11]. We assume 
that the conductivity of the background medium is 1 and that 
its permittivity is vanishing. Consider a target D — z + SB, 
where S <C 1 is the characteristic size of D, z is its location, 
and B a smooth bounded domain containing the origin. We 
assume that D is of complex admittivity k — a + iecj, with oj 
being the operating frequency in the range [1, 10] and a and 
e being respectively the conductivity and the permittivity of 
the target. It has been also shown in [I] that, in the presence 
of D, the electric field u generated by the fish is the solution 



of the following system: 

Au = f 

V-(l + (fc-l)xu)Vu = 

P =0 
dv 



u\+ -u\- 
\u(x)\ 



in Q, 

in K d \ H, 

on dtt, 
on d£l 



[1] 



= 0(|: 



Here, xd is the characteristic function of D. Fig. [T] shows 
isopotentials with and without a target with zero permittiv- 
ity but different conductivity from the surrounding medium. 
Note that if the target's admittivity depends on the frequency 
(i.e., if the permittivity is nonzero), then a phase shift in the 
electrical potential is induced. 

In a previous study [2J, we have extracted the GPTs of 
a target by multistatic measurements using arrays of sources 
and receptors. These GPTs were then arranged and compared 
to a dictionary of already known shapes. This study aims at 
adapting this method to the electro-sensing problem. 

Asymptotic formalism. The first step is to compute the GPTs 
from the measurements. In this regard, the next result will 
be useful. Except when mentioned, we will fix in this section 
the frequency to, leading to a fixed complex admittivity k. We 
will need the following notation. For a multi-index a G N d , 



let x° 



d ad 



with dj = d/dxj. Let 



G(x) be the Green function of the Laplacian in R d which sat- 
isfies AG = 5 (where <5 is the Dirac function at the origin) and 
is given by 



G(x) = 



2tt ' 



■In I x I 



1 1 

4-7T \x\ 



d = 2, 
d = 3. 



We denote the single and double layer potentials of a function 
4> G L 2 (dft) as Sn[4>] and T>n [<f>], where 

Sn[c/>](x):= [ G{x-y)^{y)do{y), x G R d , [2] 
Jan 



and 



dG 



VoMx):= I 1 f^{x~y)4>{y)da{y), x£R d \dn. [3] 
Jan ov \V) 

We also define the boundary integral operator ¥i\i on L 1 {dtt) 
by 

r fift 

K5MW:=/ ^—Jx-y^dviy), <j> G L 2 (dQ). [4] 




Fig. 1. Isopotentials without (a) and with (b) a target with o = 5 and £ = 0. 
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The operator ICq is called the Neumann-Poincare operator. 
We assume that the target is away from the fish, i.e., the dis- 
tance between the fish and the target is much larger than the 
target's characteristic size but smaller than the range of the 
electrolocation which does not exceed two fish's body lengths. 
The following theorem holds. 

Theorem 1. Let us define the function H : R d — > C by 



H(x) = p(x) + Sn 



where p is the field created by the dipole f, i.e., Ap = / inM. d . 
Then, for every integer K > 1, the following expansion holds 

K K-\a\+l 



The third line gives us 



du 






du 




dv 


+_ 


dv 


+_ 



u(x) = H{x) + 5 d 



E E 

H=i \p\=: 



xM a/3 (X,B)d ,3 G{x- 



f_ 1 )l/3| ( 5kl + l/3| 
L22 ^ daH{z) 

I 

z) + 0(S d+K ), 



[6] 

uniformly for x G dil, where z is the location of the target D 
and 

M a p(X,B)-.= [ (AJ-ZC!)- 1 ^]/^) 
Jbb dv 

is the generalized polarization tensor (of order (a,j3)) associ- 
ated with the domain B and the contrast A = (k + l)/2(fc— 1); 
see Here, K,* B is the Neumann-Poincare operator associ- 
ated with B and I is the identity operator. 

Let us make a few remarks. First, the definition of the 
GPTs still holds for complex-valued A. However, some prop- 
erties are lost by this change; thus one has to study them more 
carefully in this situation. Second, the function H, which is 
computed from the boundary measurements, still depends on 
5 but this is not important for our present study. Indeed, 
formula (J6j> could have been derived with U, the background 
solution m the absence of the target, instead of H and Gr 
- the Green function associated to Robin conditions on (9f2 - 
instead of G, but it is much easier to compute d a H(z) and 
d l3 G(x — z) once z is known. This leads us to the third remark: 
the location z is supposed to be known from the algorithm de- 
veloped in p]. Electrolocation algorithms are either based on 
a space-frequency approach in the case of multifrequency mea- 
surements or on the fish's movement if only one frequency is 
used [T| |5T] . Finally, it is worth mentioning that the knowl- 
edge of M a p{X, B) for all a, /3 € N d determines uniquely B 
and A [7] Moreover, the following scaling relation holds: 

M aP (\,SB) = 8 d - 2+M+w M a p(\,B). 

We will follow the proof of [8] Theorem 4.8]. In a first 
step, let us show the following formula. 
Lemma 2. For x G K d , 



u(x) = H{x) + S D (XI - K* D )~ 



1 


' dH 






dv 


3D- 



(*), 



[7] 



where v is the outward normal unit vector at dD and Sd and 
K* D are respectively defined by p|) and (01) with Q replaced with 
D. 

Proof. In [JQ Section 4.1.2], it is shown that 

u(x) = p(x) + Sa[iP](x) +V n [<t>](x) + S D [<f>](x), 

where the functions ip, 4> G L 2 (dQ.) and <j> G L 2 (dD) verify 
the following system 

4> — on d£l, 

2 dv J dv dv 

~(5oM)+e^(%M) + (A/-K:i)M = gjon dD. 



1 


dH 






dv 


3D. 



and the jump formulas for the single and double layer poten- 
tials Theorem 2.17] give us 

du 



dv 



and d> — —u\ 



so that, from the boundary conditions of the system |TJ, we 
obtain p + Sn [tp] + [<j>] = H and the lemma is proved! □ 

We can now prove Theorem 111 using the arguments in [8] 
pp. 72-73]. Starting with formula 17b, the proof relies on a 
Taylor expansion of H and the Greenfunction involved in the 
single layer potential. Indeed, denoting 

K 1 

Hk(x)= J2 -,d a H(z){x-z) a , 
1 — ' a! 

|a|=0 

a Taylor expansion gives us 
dH dH K 



dv 



dv 



< C5 K \dD 



1/2 



L 2 (3D) 

and from [S] Formula (4.10)], we have for any h G L 2 {dD) 
such that f gD h — 0: 

Vz G dQ, \S D (\I -tChy^mx^ < CS\dD\ 1/2 \\h\\ L2(9D) . 

Hence, using the fact that |<9D| = <5 d_1 \dB\, we obtain 



S D (XI-Ki 



< C8\dD 



1/2 



dH 


dH K 


dv 


dv 


dH 


dH K 


dv 


dv 



< CS 



d+K 



L 1 (dD) 



Plugging this inequality into Q enables us to write, for 
x G dQ, 



u(x) = H{x) + S D (XI - K*o)~ 



dHK 
dv 



(x) + 0(S d+K ) 



By a change of variables y — (y — z)/8, denoting <j> a {y') = 
(XI - £b)~> ■ Vw a ](y') for y' G dB (where v is here the 
outward normal unit vector to dB), we have (see for example 
the arguments in [5] Section 3]) 

K 1 

u(x)-H(x)= ^d a H(z)5 ]al+d - 2 

\a\=0 °' 

x f G(x-z- 8y')<t> a {y')do-(y) + 0(8 d+K ). 

J dB 

We can now conclude by injecting a Taylor expansion of the 
Green function 



G(x-Z- 5y) = Y, 
l(3|=0 

in the integrand, giving 



d^Gix - z)y^ 



u(x)-H(x) = 5 d -*Y: E L ^ L A, ^././(,) 

W=o |^|=o 



al/3 



xd p G(x-z) / y p <j> a (y)d°(y) + o(s a+K ). 

J dB 

The last term is M a p{\, B) by definition [9]; it then suffices 
to show that the terms with \a\ = or |/3| = vanish, which 
is the case because J gB <j> a — and <j> a = if |aj = 0. Thus, 
Theorem [l] is proved. □ 
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Data acquisition and reduction. Let us suppose that the fish is 
moving, and let us take a sample of S G N* different positions 
(r2 3 ) 1<s<s . This gives us 2S different functions (u 3 ) 1<s<s and 
(H s ) 1<g<s , leading us to the following data matrix 

Q:= (Q«0i< s < S ,i<r<fl == (u s (xi S) ) ~ H s (x^)) 

— — \ /1<S<S,1<?'<K 

where (i^ g 9fi s )i<r<H are the receptors of the fish being 

in the s^ 1 position. The choices of indices are motivated by 
the fact that the different positions play the role of sources. 

The goal of this subsection is to simplify this data set in 
order to extract the GPTs (precisely, the contracted GPTs as 
it will be seen below). Indeed, from pj, one has 



Moreover, Ai a \ s x and B\ a \ s , z can be expressed in the follow- 
ing way 



k K-\a\ + l cM + l/31 

|a|=l |/3|=1 



xM a0 (X,B)d p G(x { r s> -z) + 0(S d+li ). 

As in [2], we will express this formula in terms of contracted 
GPTs (CGPTs), when the dimension of the space is d — 2. 
Let us first recall the definitions of these contracted GPTs. 
For a target B with contrast ratio A, knowing the GPTs 
M a p{\, B) for all indices a and /3 such that \a\ = m and 
\j3\ = n leads us to construct the following combinations, 
called CGPTs, 

M mn = E E a ™a}M al3 , 

\ct\— m |/3|— n 

C» = E E a a hr p M af), 
\ct\— m |/3|— n 

M% n = E E b ™ a dMc>i3, 

\ct\= m |/3|— n 

ac„ = E E WpMc*, 

\a | — m | /3 1 — 71 

where the real numbers a™ and &2 1 are defined by the following 
relation 

/ . • \m \ ^ ma.- \ r jm 

(Xl + 1X2) — 2_^, a a X 0p X ■ 

|a|— m |/3|— m 

In the polar coordinates, x — r^e l8x , these coefficients also 
verify the following characterization 

\ mam rt l \ r 1 m 8 m n 

> a a x = r x cosmo x and y, °/3 x = r x smmo x . 

\a\=m |/3|— m 

This enables us to show 2, Appendix A. 2] that 



1 d a G{x) 



2tt \a 



1 / | a | cos \a\ 6 x ,\ a \sxa\a\6 x 
1 a'„ — ho' 1 



.1" 



[10] 

From the definition of H, and with the help of the previous 
formula, one can prove the following lemma. 
Lemma 3. Let the source f be a dipole of moment p s placed 
at z s : 

Ps (i) = Ps -VG(i-4 [11] 
Then, for any a G N 2 , there exist two real numbers A\ a \ tS ,z 



1 d a H s (z)=A }aUs , z a^+B l , 



A-T, 



(-iy 

2tt 



4>m+l{z — Zs) 
1pm+l(z - Zs) 



an dv 

iff 4>m+i{y - z) \ du 



1 

2irm 



2tt 



/ 

J an 



(y)<t>m{y - z)do{y), 
V 

Tpm + l(y 



v v (y)do(y), 
OV j 



(-ir 

2tt 



ipm + \{z — Z s ) 
— <t>m+l{z — Z s 



1 f du s 
2™ J an dv 



(y)i>m(V ~ z)do{y) 



2tt 



where the functions 4>m and tpm are defined for x G R 2 , 
x = (r^ftj) in the polar coordinates, by 

a ( \ — cosm ®x i ( \_ smm9 x 

<f>m(X) — — , 1pm(X)— — • 



Proof. Let us fix a G N 2 and define m — \a\. Let us recall the 
definition of H, given in |5| 



H a (x) =p s (x) +Sn B 



du a 










dv 


+_ 


dv 


+_ 



where Ap s = / s in R . From (111 it follows that 



d a p B (x) = p s ■ Vd a G(x — z s 



Hence, ( 10 1 yields 
(_!)!«■ 



a Ps(z) = a a 



2irm 



Ps • V0 m (2 - Z s ) 



2nm 



Moreover, we have 

In the same manner, from 
du. 



Ps • Vljjm(z — Zs) 



Vip m = —m 



Tpm + 1 
— <j>m+l 



Sn s 



we can deduce 



-,d a Sn s 
a). 



(x) 



du s (y) 



on. 



dv 



G{y - x)do(y), 



a m ( [ — 



du s 




dv 


+_ 


du s 




dv 


( 

+ 




1 



du a 
2irm dv 



\a\aa, 



, > cos m6y~ z 

(y) do(y) 

+ y-z 



+ ° a 



1 du s i sSmm6y- z . . 

(v) „m da (y) 



on, 



2-um d, 



Combining those two equations leads us to the desired result. 
□ 
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From (|9j> , the data matrix is then expressed as follows 



K+l 



Qrs — ^ ^ ^^la^s^^L ' "T" -^|a|,s,2^L 'J 
|a| + |/3|=l 

■ bl p 1 sin 

x M Q/3 (A,<5B)- 



: /3 

101 



2tt LS r" 



+ 0{S K+2 ) 

K+l 

m+n = l"" " 



cos n6 
sin ?!# 



(«) 



27T7ir" 



GW, 



Thus, defining the following matrices 



[12] 



( Mn M12 

M21 



M 1K \ 




E = (E rs ) 1 



<r<R,l<s<S ' 



\ Mai ... / 

the problem is to recover the matrix M knowing the matrix 
Q = £(M) +E, 



where £ is the linear operator defined by ( 12 1. 

Therefore, the CGPTs of the target JJ can be recon- 
structed as the least-squares solution of the above linear sys- 
tem, i.e., 



argmm | 

M±ker(£) 



£ 



[131 



where ker(£) denotes the kernel of £ and || • \\f denotes the 
Frobenius norm of matrices [2] . 

Let us remark that, in the case of multifrequency measure- 
ments (u\, . . . ,u>f), we can reconstruct [ M (/) ) from 

V ) l<f<F 

(Q^-M analogously. 

V / l<f<F 

It is worth mentioning that the location of the target de- 
tected by the fish may be different from the true one. More- 
over, the target may be rotated and hence, the reconstructed 
CGPTs correspond to a translated, scaled, and rotated target 
B. In order to recognize the shape B, it is therefore funda- 
mental for the recognition procedure to have size invariance, 
rotational invariance, and translational invariance. This could 
be related to the behavioral experiments that have shown 
that weakly electric fish categorize targets according to their 
shapes but not according to sizes, locations, or orientations 
El- 



Recognition and classification 

Depending on whether we consider multifrequency measure- 
ments or not, we will not identify the CGPTs in the same 
way. 

Fixed frequency setting: shape descriptor based classifica- 
tion. In [2], an algorithm based on shape descriptors was de- 
veloped for the recognition of a target in a more classical elec- 
trical sensing setup, (i.e., multiple sources/receptors placed on 
the surface of a disk). In this paper, we apply this algorithm 
in the context of electro-sensing. 



We recall here the concept and properties of shape de- 
scriptors in two dimensions [2]. For a double index mn, with 
m, n = 1, 2, . . ., we introduce the following complex combina- 
tions, N« = (jV^ n ) m ,n,N^ = (Ad 2 2) m ,n, of CGPTs: 



A/"™ (A, D) = (M- n - AO) + i{MZ n + M^ n ), 
N£l(\, D) = (AC„ + AC«) + i(M^ n - AC„). 



[141 



Let 



N*${D) 



2N${D) 



, T- u D = {x~u,xeD}. 



We define the following quantities which are translation in- 
variant: 

J W {D) = N (1) (r_„£>) = C-"N (1) (D)(C- U )', [15] 

J (2) (D) = N (2) (r_„D) = C^N (2) (D)(C~ U )', [16] 

with t being the transpose and the matrix C~" being a lower 
triangular matrix with the m, n-th entry given by 



1 \ r 
(-«) 



From^'p) = (j£l(D)) m ,n and J^(D) = {j£l(D)) m , n , 
we define, for any indices m, n, the scaling invariant quanti- 
ties: 



\J mn (-^) 



1/2 ' 



1,2. 



[17] 



Finally, we introduce the CGPT-based shape descriptors 

jw _ ( T m ) anf1 T (2) _ f T (2) 



l£l(D) = \S<»(D)\, l£l(D) = \S£l(D)\, 



[181 



where I ■ I denotes the modulus of a complex number. Con- 
structed in this way, and X^ are invariant under trans- 
lation, rotation, and scaling. Only shape descriptors of order 
2, i.e., for m,n = 1,2 will be used in the sequel. Shape de- 
scriptors in three dimensions were derived in [4]. 

Multifrequency setting: Spectral induced polarization based 
classification. When multiple frequencies are involved, we can 
use the shape descriptors Imi(fl) and Xmn(D) at frequencies 
cji, . . . ,ujp to enhance the stability of the classification. How- 
ever, as it will be shown later, this does not yield a very stable 
classification procedure. 

Here we rather focus on the first-order polarization tensor 
(PT), that is, the 2x2 complex matrix M^\D) associated 
with the target D and frequency /: 



M U \D) 



a + 1 + iujfe 
an V 2 ^ -1 + iw/e)' 



fCo) [u]ydcr(y), 



for / = 1, . . . ,F. We will show that they are sufficient to 
identify efficiently the targets. Note that it is not possible 
to compute the shape descriptors Xmn(D) and Xmn(D) based 
only on first-order PT, because they require at least second- 
order polarization tensors. This limits the use of shape de- 
scriptors in the limited- view case where the reconstruction of 
higher-order GPTs is not accurate [3]. 

Here we use the spectral content of the first-order PTs for 
recognition. We have the following properties [5]. 
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Proposition 4. For any scaling S > 0, rotation angle 6 £ 
and translation vector z £ ffi 2 , let us denote 



D = 2 + (5i? 9 B := {x = z + <5i? s 



where 



cos t/ — sin ( 
sin 6 cos 



% := 

is the rotation matrix of angle 6. Then, 



M U) {D) = S 2 R e M u) (B)Rj. 



[191 



Hence, if we denote by t£ (D) and t 2 (D) the singular val- 
ues of M^ '(D), we obtain 



Vie {1,2}, rY\D) = S\. 



This gives an idea for two algorithms: 



The first one, matching the singular values of all the first- 
order PT ( My' J , would be dependent of the char- 

V / x<f<F 
acteristic scale 8 of the targets in the dictionary; 
The second one, independent of the scale of the target, 
would match the following quantities 



for i = 1,2 and / = 1, 



^ = 



;F 



JF) ' 



[201 



Some comments are in order. First, the reason why we con- 
sider the first one, even if it is scale-dependent, is because it 
is far more stable. Also, in some biological experiments, two 
targets of different scales are considered as different [35]. A 
question raised was then: how is it possible to discriminate 
between a nearby small target and an extended one situated 
far away? With the second algorit hm , we have an answer. 
The last remark concerns equation (20 1. We could have also 



considered other scale-dependent ratios, such as 



Jf) 



J 1 ) 



Jf) 



but since rj F ' happens to be the largest one (the frequen- 
cies are sorted in increasing order), it is more stable to con- 
sider n20n. It is worth mentioning that if there exists an integer 
p > 2 such that R 2w / p D = D, then M^\D) is proportional 
to identity. 

Background field elimination. We can also improve stability 
of reconstruction by eliminating the background field. Let us 
denote by U(x) the background electric field {i.e., the solution 
of |l| with k = 1). In [T| Proposition 2], we have proved the 
following formula: 



Vv. 



dv 



dU 
dv 



VU(z) T M (f) 



(£)V, 



8G 

dv x 



[21] 

where Uf is u associated with the / th frequency, M^\D) is 
the first-order polarization tensor at the f th frequency, and 
Vn is the (real-valued) postprocessing operator given by 

-d ■ l T r* t dVn 
Vn.= 2 J -^-^' 



with Vn and /Co being defined by (Bl and Q; see [T]. Hence, if 
the emitted signal U is real-valued7~then taking the imaginary 
part leads us to 



dv 



W( 



"QmM u) (D)X7 z 



dG 
dv x 



[22] 

Note that in biological sciences, the restriction on U to be real 
is justified since the per mitt ivities of water and the fish are 
negligible [I] . Now, from ( 22 1, we can extr act Qm M ^ (D) by 
solving a least-squares problem similar to ( 13 1. Then, we have 



the singular values of the imaginary part of M U) (D), which 
would be sufficient for shape recognition and classification. 
The goal of th is procedure is to get rid off the computation of 
dU/dv in \21\ , which is supposed to be performed numerically 
in real- world applications, thus subject to errors. Note that 
the postprocessing operator Vn makes the data independent 
of the shape of the fish's body. 





Fig. 2. The 8 elements of the dictionary. The dotted lines indicate a target with 
different electrical parameters. 





~ . . ./ 












Fig. 3. Two different kinds of experiences involve (on the top) a twisted-ellipse 
shape or (on the bottom) ellipse shape. The real part of the electric field is plotted, 
for 3 (out of 20) positions that the fish takes around the target (placed at the origin). 




Fig. 4. A fish having the shape of a twisted ellipse with a smaller aperture than 
in Fig.fJ 



6 | www.pnas.org 



Footline Author 



Because of the following relation which follows from ( 19 1 
M U) (D) = 0(S 2 M (f) (B)), 



frequencies, equally distributed from ojq := 1 to Fojo- Again, 
we refer to [T] for nondimensionalization of the underlying 
model equations with proper quantities and realistic values. 



taking the imaginary part would lead us to only compute 
X7U (z) in ( |22[ ) and hence, the error made here would be mod- 
ulated by a factor of order S 2 . 

Numerical illustrations 

In this section, we illustrate the performance of the algorithms 
developed in the previous section. We use the CGPTs ob- 
tained in order to classify the targets. We present an example 
with fixed frequency, and another with multifrequency mea- 
surements. As it will be seen, the latter does not lead us to a 
significantly more stable classification in the presence of noise 
or for limited aperture. The errors in the reconstruction of the 
high-order polarization tensors due to measurement noise or 
the limited-view aspect deteriorate the stability of the pro- 
posed algorithm. However, when, at multiple frequencies, 
only the first-order polarization tensor is used, we arrive at 
a very robust and efficient classification procedure. 

For the sake of clarity, and due to the large numbers of 
computations performed, the results are presented in the ap- 
pendix. 

Setup and methods. We describe the dictionary as well as the 
measurement systems. We consider two different shapes for 
the fish: ellipses and twisted ellipses. Note that this vari- 
ety of shapes exists in nature. On the one hand, twisted el- 
lipses would represent electric eels (Electrophorus electricus), 
whereas on the other hand straight ellipses would look like 
Apteronotids [35]. This simplified representation shows that 
the principle of our algorithms can be generalized to any kind 
of fish's shape (hence modeling, for example, electro-sensing 
for Mormyrids as well). It also enhances the fact that, for 
bio-inspired engineering applications, the shape of the robot 
is not determining. Moreover, as we will see later, our simpli- 
fied representation is a good model to tackle aperture issues. 

Dictionary 

The dictionary T> that we consider is composed by 8 differ- 
ent targets: a disk, an ellipse, the letter 'A', the letter 'E', a 
rectangle, a square, a triangle, and an ellipse with different 
electrical parameters (see Fig. [2). Indeed, all the other tar- 
gets have conductivity a — 2 and permittivity e = 1 whereas 
the last one has conductivity a = 5 and permittivity e = 2. 
Except when mentioned, the characteristic size of the target 
will not matter, and will be fixed to be of order 1. 

Measurements 

In each numerical experiment, one target of the dictionary is 
placed at the origin, while the fish swims around it. As it 
has been mentioned, we consider two different shapes for the 
fish's body: ellipses and twisted ellipses. The measured data 
is built taking 20 positions of the fish around the target (see 
Fig. §. 

In Fig. HI a smaller aperture than the one in Fig. p] is 
considered. 

The typical size of the target is 8 = 0.3 while the fish is 
turning around a disk of radius R = 1; the twisted ellipse's 
semi-axes are a = 1.8 and b = 0.2 while the straight ellipse's 
semi axis are a = 1 and b = 0.2. The effective thickness of 
the skin is set at £ = 0. The fish has 2 7 receptors uniformly 
distributed on its skin, and the electric organ emits F = 10 
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Classification 

The recognition process is as follows. When measurements are 
acquired, we perform least-square reconstruction of the (first- 
or second-order) CGPT of the targets. From this CGPT, we 
compute quantities of interest q (i.e. Shape Descriptors or 
singular values or ratios of singular values). Then, for each 
element n in the dictionary T>, we compute \\q — q n \\, where q n 
is the - pre-computed - quantity of interest for the nth shape. 
This leads us to charts such as the ones presented in Figs. [5] 
§ and[7| 

Framework for algorithms of multifrequency classifica- 
tion: 

i: Input: the quantities of interest [q^') calcu- 

V / l<f<F 

lated from the measurement of an unknown shape D; 
2: for B n G T> do 

3: e n <- E!</< F \\q{B n ) U) -q U) f where (<?(£„)<«) ^ 
is the same type of quantities of interest of the shape 

4: n •<— n + 1; 
5: end for 

6: Output: the true dictionary element n* argmin n e„. 
Stability analysis 

First, let us explain what kind of noise is considered. We 
will add a random matrix (with Gaussian entries) to the data 
matrix Q defined in (JsJ) . More precisely, we will consider 

Q := Q + eW, 

where W is a S x R matrix whose coefficients follow a Gaussian 
distribution with mean and variance 1. The real number e 
is the strength of the noise, and will be given in percentage 
of the fluctuations of Q, (i.e., max Sir Q sr — min S;r Q Br )- The 
recognition procedure remains the same. 

Stability analysis was then carried out empirically: for 
each noise level, we performed JV s t a bii independent recognition 
process (with iVstabu being precised for each experiment), and 
computed the ratio of good detection. The computation ends 
when we reach the threshold of 12.5% probability of detection 
that corresponds to a uniform random choice of the object. 
This gives us Figs. [8lto[T5| 

Results and discussion. In this subsection, we compare the re- 
spective stability of our algorithms. Due to the large number 
of situations and computations, only significant figures were 
plotted, giving the following classification of recognition pro- 
cesses, according to their range of application. 

Fixed frequency setting: shape descriptors 

If only one frequency is accessible for the measurements, then 
the only algorithm possible is classification based on shape de- 
scriptors. Indeed, first-order PTs are not enough to discrim- 
inate between objects. However, the use of shape descriptors 
is limited to the twisted-ellipse case with nearly full aperture 
(see Figs. [6] and [7] where some targets are not recognized, 
such as the ellipse m Fig. [fj] or the disk in Fig. [?J. Moreover, 
its stability with respect to measurement noise is quite poor 
(see Fig. [8b. 
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Multifrequency setting: spectral content of PTs 

In the case of multifrequency measurements, shape descrip- 
tors do not increase their stability enough compared to sin- 
gular values (see Fig. |9J. Henc e, it i s better to use singular 
values of the PTs (see 



lg.S. 



10 to 13 1. One can see that: 



considering all singular values (Figs, 
more sta ble tha n considering ratios 
(Figs. 



10| and 111 is much 
"of singular values 



12jand 13|; 

the aperture does not change very much the stability. 

In this regard, the most stable algorithm is to recognize all 
singular values when the fish is a twisted ellipse (Fig 



10) 



leading us to a probability of detection superior to 90% with 
noise level of 125%. This is a huge gap when compared to 
the recognition process with shape descriptors, allowing only 
a few percents of noise. Note that the noise level is computed 
with respect to the perturbation in the measurements Q given 
by |8J, which is of order of the target volume, see ([91. Hence, 
a noise level of 125% remains small compared to the actual 
transdermal potential u. 

Background field elimination 

We can see in Figs. [14] and [15] that taking the imaginary part 
of the measurements in order to avoid the computation of the 
background field does not significantly decrease the stability 
of the reconstruction based on spectral content. Since the 
reconstruction of CGPTs is very fast, the background field 
elimination technique would yield to real-time imaging. 



Concluding remarks 

In this paper, we have successfully exhibited the physical 
mechanism underlying shape recognition and classification in 
active electrolocation. We have shown that extracting gener- 
alized polarization tensors from the data and comparing in- 
variants with those of learned elements in a dictionary yields 
a classification procedure with a good performance in the full- 
view case and with moderate measurement noise level. How- 
ever, this shape descriptor based classification is instable in 
the limited-view case and for higher noise level. When mea- 
surements at multiple frequencies are used, the stability of 
our classification approach is significantly improved by using 
phase shifts and keeping only the first-order polarization ten- 
sor. The resulting spectral induced polarization based classi- 
fication is very robust. 

Our results open the door for the application of the ex- 
tended Kalman filter developed in [3] to show the feasibility 
of a tracking of both location and orientation of a target from 
perturbations of the electric field on the skin surface of the 
fish. It also remains to understand to what extent the spec- 
tral induced polarization approach could help us retrieve the 
electric parameters of the target or locate and recognize mul- 
tiple targets. 
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Fig. 7. Results of the matching with the dictionary for a twisted ellipse shaped 



fish with smaller aperture (see Fig. |4*|. 



36. G. Von der Emde, S. Schwarz, L. Gomez, R. Budelli, and K. Grant. Electric fish 
measure distance in the dark. Science, 260:1617-1623, 1993. 



Appendix 

In this appendix, we numerically illustrate the main find- 
ings in this paper and show the potential of electro-sensing 
for shape recognition and classification. 




Shape Label 

Fig. 5. Results of the matching with the dictionary (differences between shape 
descriptors), when the fish is a twisted ellipse. The red bar represents the bar of the 
target being identified. In the ^-coordinates, 1 stands for the disk, 2 for the ellipse, 3 
for the letter A, 4 for the letter E, 5 for the square, 6 for the rectangle, 7 for the trian- 
gle and 8 for the ellipse with different electrical parameters. In the y-coordinates, the 
distance between the shape descriptor of the target - computed from measurements - 
and the shape descriptors of the dictionary. 
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Fig. 8. Stability of classification based on Shape Descriptors, when the fish is 
a twisted ellipse. Here, only one frequency (the smallest one) is considered. The 
threshold of 12.5% that corresponds to a randomly chosen target is represented in 
red dotted line. Here, N stah u = 10 5 . 
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Fig. 9. Stability of classification based on multifrequency Shape Descriptors, when 
the fish is a twisted ellipse. Here, iVst a bil — 5.10^. 
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Fig. 6. Same results as in 5 with an ellipse-shaped fish. 
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Fig. 10. Stability of classification based on differences between all singular values 
of PTs, when the fish is a twisted ellipse. The characteristic size of the target is 
supposed to be known. Here, iV s t a kii = 5-10 . 



Fig. 13. Stability of classification based on differences between ratios of singular 
values, when the fish is a straight ellipse. Here, A^ s ^ a bil — 5.10^. 
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Fig. 11. Stability of classification based on differences between all singular values 
of PTs, when the fish is a straight ellipse. The characteristic size of the target is 
supposed to be known. Here, iV s t a bil = 5.10^. 



Fig. 14. Classification with imaginary part of the PT, when the fish is a straight 
ellipse. All singular values are considered to discriminate between the targets. Here, 
iVstabil = 10 5 . 
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Fig. 12. Stability of classification based on differences between ratios of singular 
values, when the fish is a twisted ellipse. Here, -/V s t a bil — 5.10 . 



• A 
t E 

- Souare 
Recangie 

- Tnangle 

• DiHe'eni eii'pse 
10% 



Strength ot m 



Fig. 15. Classification with imaginary part of the PT, when the fish is a straight 
ellipse. Ratios of singular values are considered to discriminate between the targets. 
Here, 7V stabil = 10 5 . 
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